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Abstract. We study some function-theoretic properties on complete smooth 
metric measure space {M,g,e~"dv) with its Bakry-Emery Ricci curvature 
bounded below. Assuming the bound of function /, we derive a Moser's para- 
bolic Harnack inequality of the /-heat equation, which leads to upper and lower 
Gaussian bounds on the /-heat kernel. We also prove various L p -Liouville the- 
orems in terms of the lower bound of Bakry-Emery Ricci curvature and the 
bound of function /, which generalize the classical Ricci curvature case and 
the ^-Bakry-Emery Ricci curvature case. 



1. Introduction and main results 

1.1. Background. Let (M,g) be an ^-dimensional complete noncompact Riem- 
mannian manifold and / be a smooth function on M. We define a symmetric 
diffusion operator Af (or /-Laplacian), which is given by 

A f := A-V/-V, 

where A and V are the Laplacian and covariant derivative of the metric g. f- 
Laplacian Af is the infinitesimal generator of the Dirichlet form 



£((/> 1 ,<h) = / (V<£i,V<fe)d//, V<^ 2 £ C °°(Af), 

JM 

where \i is an invariant measure of Ay given by d/i = e~'dv, and where dv is the 
volume element induced by the metric g. It is well-known that Ay is self-adjoint 
with respect to the weighted measure d/i. The triple (M,g,e~-'dv) is customarily 
called a complete smooth metric measure space. On this measure space, we often 
consider the /-heat equation 

(d t -A f )u = 

instead of the classical heat equation. If the function u is independent of time t, 
then u is the /-harmonic function. In this paper, we denote H(x, y, t) be the /-heat 
kernel. That is, for each x € M, H(x, y, t) — u(y, t) is the minimal solution of the 
/-heat equation with u(0,y) = S x (y). Equivalently, H(x,y,t) is the kernel of the 
semigroup P t = e tAf associated to the Dirichlet from £(<f), (f). 

On the metric measure space (M, g, e~^dv), following Bakry and Emery [T] and 
[2] (see also [22] and [24]), we define the Bakry-Emery Ricci curvature as follows: 

Ricf := Ric + Hess(f), 
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where Ric denotes the Ricci curvature of the manifold and Hess denotes the Hessian 
with respect to the manifold metric. The Bakry-Emery Ricci curvature is a natural 
extension of Ricci curvature. If / is constant, Ricf returns to the Ricci curvature 
Ric. The Bakry-Emery Ricci curvature has been extensively studied because it 
often shares similar properties of the ordinary Ricci curvature. For example, there 
exists an interesting Bochner type identity 

A/|Vu| 2 = 2| J ffess(u)| 2 + 2(Vu, VA/u) + 2Ric s (Vu, Vu). 

This identity is parallel to the Ricci curvature case, which plays an important role in 
studying comparison theorems (see [3U])- For more extended results, the interested 
reader can consult [3], g], [SJ, 03, [23], [2J], [25], [35], [36], [38] and gD]. 

Also, the Bakry-Emery Ricci curvature has become an important object of study 
in the geometric analysis, in large part due to the so-called gradient Ricci soliton. 
Recall that a complete manifold (M, g) is a gradient Ricci soliton if 

Ricf = pg 

holds for some function / and constant p. The soliton is called expanding, steady 
and shrinking, accordingly, if p < 0, p = and p > 0. Ricci soliton of manifolds 
possess many interesting geometric and topological properties. See, for example, 
[5J, [7J and [2U] for nice explanations on this subject. 

Recently, there has been renewed interest in the Bakry-Emery Ricci curvature. 
For example, G. Catino, C. Mantegazza, L. Mazzieri and M. Rimoldi [9], P. Peterson 
and W. Wylie [30], and S. Pigola, M. Rimoldi and A.G Setti [31] established various 
Liouville-type or rigid theorems on this curvature. Prior to their works, X.-D. Li |22] 
studied the L 1 -Liouville theorem if the A-Bakry-Emery Ricci curvature is bounded 
below by a negative quadratic polynomial of the distance function, which is an 
extension of the classical Z^-Liouville theorem on Ricci curvatures, proved by Peter 
Li [18]. Here the A-Bakry-Emery Ricci curvature is defined by 

Ricf :=Ri c + Hess(f)-^j^-, 

where A is a positive constant. However, as X.-D. Li said in Subsection 8.6 of 
[22], we cannot prove the Z^-Liouville theorem if we only keep the lower bound 
of the Bakry-Emery Ricci curvature to be a negative quadratic polynomial of the 
distance function. Because we can not obtain the Li-Yau type parabolic Harnack 
inequality under only this curvature assumption. Now it is natural to pose the 
following problem: What is the optimal geometric and analytic condition on the 
metric measure space such that the Li-Yau parabolic Harnack inequality holds? 

In recent papers [27] [25], Munteanu and Wang partially answered the above 
question. In particular, they derived gradient estimates and Liouville properties for 
positive /-harmonic functions under some growth assumption on /. Their theorems 
take the form of Yau's classical result on the positive /-harmonic functions, but 
their proof is new and quite different in spirit from Yau's direct Bochner formula 
method [52] . Their proof approach essentially relies on the well-known De Giorgi- 
Nash-Moser theory, which motivates our proof of Theorem 11.11 in Subsection 11.21 

1.2. Main results. The purpose of this paper is to study some geometric inequal- 
ities for the /-heat equation and various L p -Liouville theorems for the /-harmonic 
function on complete smooth metric measure space (M,g,e~^dv) with its Bakry- 
Emery Ricci curvature bounded below, which is divided into two parts. 
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In the first part, borrowing the idea of Munteanu and Wang in [27l EH] , we will 
study some geometric inequalities, such as Harnack inequalities, Holder continuity 
estimates, /-heat kernel estimates on complete smooth metric measure spaces. We 
first present a parabolic Harnack inequality of the /-heat equation on complete 
smooth metric measure spaces. 

Theorem 1.1. Let (M, g, e~f dv) be an n- dimensional complete noncompact smooth 
metric measure space. If Ricf > — (n — 1)K and \f\(x) < A for some nonnega- 
tive constants K and A, then there exist a constant c(n, A) such that, for any ball 
B {r), o G M , < r < R < oo and for any smooth positive solution u of the f -heat 
equation in the cylinder Q — B (r) x (s — r 2 , s), 

SUp{ U } < e ^,A){l+K^) . inf{u}; 

Q- Q+ 
where Q_ := B Q {\r) x (s - |r 2 , s - \r 2 ) and Q + := B D (\r) x (s - \r 2 7 s). 

The sketch proof of Theorem 11.11 will be given in Section [3] Its proof follows 
the Moser iteration technique [2B], which involves a local Sobolev inequality on a 
smooth metric measure space. Munteanu and Wang [28] used similar technique to 
derive an elliptic Harnack inequality for /-harmonic function. When the metric 
measure space is a Riemannian manifold, that is function / is constant, this case 
was obtained independently by L. Saloff-Coste [32] and A. Grigor'yan [14] . 

A standard consequence of Theorem 11.11 is a strong Liouville theorem for any 
/-harmonic function (see Corollary 13.21 in Section [3]). Moreover, Theorem ll.il also 
implies two-sided /-heat kernel bounds with respect to smooth metric measure 
spaces. This result is essential parallel to the case of heat equation on Riemannian 
manifolds, obtained by L. Saloff-Coste |34) . 

Theorem 1.2. Let (M,g,e~ f dv) be an n-dimensional complete noncompact smooth 
metric measure space. If Ricf > — (n — 1)K and \f\(x) < A on the ball B a (2R) 
for some nonnegative constants K and A, then there exist positive constants Ci, 
i = 5, 6, 7, 8, depending only on n and A such that 

e- c «( 1+xt ) / d 2 (x,y)\ TTI . e ^+Kt) / #{x,y) 

= — exp -c B < H(x,y,t) < ^- exp -c 7 ^ y> 

V f (B x (Vi))) P \ t J- { > U ' >- V f (B x (Vi)) P \ t 

for any x, y G B a (R/2) and < t < R 2 /4, where Vf(B x (y/t)) denotes the f -volume 
of the ball B x (y/i) with respect to e~^dv. 



Remark 1.3. Theorem 11.21 gives a detailed description of the coefficients of two- 
sided Gaussian bounds on the /-heat kernel. This estimate has an important ap- 
plication in the proof of the following Theorem 11.61 



We point out that the proof of Thcorcm ll.2l is different from the classical Li-Yau 
trick [18] . In [18] , two-sided Gaussian bounds on the heat kernel are obtained by the 
Li-Yau gradient estimate. However, in our case it seems to be impossible to utilize 
Li-Yau gradient estimate method directly to derive upper and lower bounds on the 
/-heat kernel on complete smooth metric measure spaces. In our setting, two-sided 
Gaussian bounds on /-heat kernel strongly rely on the Moser's parabolic Harnack 
inequality and the integral estimate of the /-heat kernel due to E. B. Davies |11) . 
which are similar to the arguments of L. Saloff-Coste 32, 33, 34 and A. Grigor'yan 
[14] . Please see the Section |4] for a detailed discussion. 
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In the second part of this paper, we will investigate various L p -Liouville the- 
orems for /-harmonic functions on complete noncompact metric measure space 
(M,g,e~'dv) under different assumptions on Ricf and /. 

We first start with a L p -Liouville theorem for positive /-subharmonic functions 
when I < p < oo. 

Theorem 1.4. Let (M,g,e~fdv) be an n- dimensional complete smooth metric 
measure space. For any 1 < p < oo, there does not exist any nonconstant, nonneg- 
ative, L p {ii)-integrable f-subharmonic function. This constant must be zero if M 
has infinite f -volume. 

The proof of this result is almost the same as the arguments of the classical case 
[43] , which will be given in Subsection 16.11 Note that this theorem holds without 
any curvature condition. 

We also deal with the L p -Liouville theorem in case of < p < 1. This case is 
parallel to the case of subharmonic functions on a manifold with Ricci curvature, 
discussed by Li and Schoen [20] . See Subsection 16.21 for the detailed discussion. 

Theorem 1.5. Let (M, g, e~-> dv) be an n-dimensional complete noncompact smooth 
metric measure space. If 

Ricj > S(n)(l + r{x)y 2 , 

where r(x) is the distance to a fixed point o € M , and f is bounded, then any 
nonnegative LP -integrable (0 < p < 1 ) f-subharmonic function must be identically 
constant. Moreover, this constant must be zero if M has infinite f -volume. 

Finally, motivated by the P. Li's work [18] and X.-D. Li's generalization [22], we 
obtain a new T^-Liouville theorem on the smooth metric measure space. 

Theorem 1.6. Let (M, g, e^^ dv) be an n-dimensional complete noncompact smooth 
metric measure space. Assume that 

Ricf > -C(l + r 2 {x)) 

for some constant C > 0, and f is bounded. Then any nonnegative L 1 ([i) -integrable 
f-subharmonic function must be identically constant. Moreover, this constant must 
be zero if M has infinite f -volume. 

Theorem ll.6l partiallv answers a question posed by X.-D. Li (see Subsection 8.6 in 
[2"2"] ) . Its proof is similar to the arguments of [TB] , where a critical step is the usage 
of the upper Gaussian bound on the /-heat kernel (Theorem 11.21) . The detailed 
discussion shall be carried in Subsection 16.31 

Remark 1.7. We remark that the absolute value of a /-harmonic function is a 
nonnegative /-subharmonic. Therefore we can conclude that a complete metric 
measure space does not admit any nonconstant L p (/u)-integrable /-harmonic func- 
tion under the same hypotheses of Theorems 11.41 11.51 and 11.61 respectively. 

Remark 1.8. As many recent authors said in JT2], [37] and [H], if the condition 
on / bounded is replaced by |V/| bounded, then similar results of Theorems 11.51 
and 11.61 can be immediately obtained by modifying the arguments of |22| . Indeed, 
the conditions Ric/ > — (n — \)K and |V/| < a imply that 
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In the end of this section, we would like to say that our basic idea of proving 
L p -Liouville theorems was inspired by the classical work [43] , [20] , [18] and [22] , but 
we need to do a lot of work in order to develop this basic idea in our situation. For 
example, an upper bound of the /-heat kernel is quite essential in proving the L 1 - 
Liouville theorem for /-harmonic functions. However, to the author's knowledge, 
it seems to be impossible to directly apply Li-Yau gradient estimate to derive the 
upper bound on the /-heat kernel under our curvature assumptions. Fortunately, 
we can apply Moser's parabolic Harnack inequality to derive it. 

The rest of this paper is organized as follows. In Section [5] we recall the Neu- 
mann Poincare inequality and the local Sobolev inequality on complete smooth 
metric measure spaces. After that, we establish a Moser's version of parabolic Har- 
nack inequality. In Section [3] using the parabolic Harnack inequality, we obtain a 
quantitative Holder continuity estimate for a solution to the /-heat equation, which 
implies a strong Liouville theorem. In Section 0] we prove the two-sided Gaussian 
bounds on the /-heat kernel on complete smooth metric measure spaces. In Section 
[5] we prove a mean value inequality on complete smooth metric measure spaces, 
which is similar to the case of harmonic functions on a manifold, obtained by Li 
and Schoen [20] . In Section HI we derive L p -Liouville theorems of /-Laplacian on 
complete smooth metric measure spaces by following the ideas in [IS] , [20] and [J3] . 



2. Poincare, Sobolev and Harnack inequalities 

Let Af — A — V/ • V be an /-Laplacian with an invariant measure d/i = e~^dv 
on a complete Riemannian manifold. For a set 51, we will denote by V(Q) the 
volume of with respect to the usual volume form dv, and Vf(Q) the /-volume of fl 
with respect to e~'dv. Throughout of this section, we will assume 

Ric f > -{n - l)K and \f\(x) < A 

for some nonnegative constants K and A. Under these assumptions, we have the 
Laplacian and volume comparison results with respect to the smooth metric mea- 
sure space {M,g,e~* dv). 

Lemma 2.1 (Wei and Wylie [39]). Let (M, g, e _/ cfo) be an n- dimensional complete 
noncompact smooth metric measure space. If Ricf > — (n— 1)K and \f\[x) < A 
for some nonnegative constants K and A, then 

A f r(x,y) < (n - 1 + <L4) VK coth \[Kr. 

for any < r < R. Hence along any minimizing geodesic starting from x G B {R) 
we have 

(9U Vf(B x (r 2 )) V£ +iA (r 2 ) 

( ' ) V f (B x (n)) ~ ^ +4A (n) 

for any < r% < r2 < R. Here V^ +4A (r) be the volume of the radius r-ball in the 
model space M^ +4A , the simply connected model space of dimension n + A A with 
constant curvature K . 

From (|2.ip . we easily deduce that 
(2.2) V f {B x (2r)) < 2 n + 4A e c ^ A)VWr ■ V f {B x {r)) 
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for any < r < R. This inequality possesses the volume doubling property, which 
plays a crucial role in our paper. Recall that a manifold admits a volume doubling 
property if for any fixed < R < oo, there exists a constant C such that 

V f (B x (2r))<C-V f (B x (r)) 

for any < r < R and x G M. By Lemma |2~T| following the L. Saloff-Coste's proof 
(see Theorem 5.6.5 in [34]) step by step, we can easily get a Neumann Poincare 
inequality with respect to smooth metric measure spaces. 

Lemma 2.2. Let (M,g,e~^dv) be an n- dimensional complete noncompact smooth 
metric measure space. If Ricj > — (n — 1)K and |/|(x) < A for some nonnegative 
constants K and A, then for any x € B a (R) 

(2.3) f \cp-cp Bo{r) \ 2 e-fdv < e Cl(1+VWr) ■ r 2 f \V V \ 2 e- f dv 

JB (r) JB„(r) 

for all < r < R and ip £ C°°(B (r)), where <fB (r) '■— Vf (-^°( r )) J B r r \ <pe~^dv. 
The constant c\ depends only on the dimension n and A. 

Remark 2.3. Munteanu and Wang [27] proved the above Poincare inequality under 
a weaker hypothesis on the oscillation of / on unit balls based on the arguments in 
[5]. That is, (|2.3[) holds on the unit balls if Ricf > —{n—l)K and sup ygSx ( 1 ) \ f(y) — 
f(x)\ < a for any x € M. But in our case, the Poincare inequality can hold on 
any radius of balls due to a stronger assumption on /, which is a crucial step on 
the proof of L 1 -Liouville theorem. Because in the course of proof of L 1 -Liouville 
result, we need to let the radius of balls tend to infinity. Also note that when / is 
constant, (|2.3|) was confirmed by L. Saloff-Coste (see (6) in [33] or Theorem 5.6.5 
in [31]). 

Combining Lemma 12.11 Lemma 12.21 and the argument in [32] , we have a local 
Sobolev inequality, which is one of the key technical points needed to apply Moser's 
iterative technique to derive parabolic Harnack inequalities for the /-heat equation. 

Lemma 2.4. Let (M,g,e~^dv) be an n-dimensional complete noncompact smooth 
metric measure space. If Ricf > — (n — 1)K and |/|(x) < A for some nonnegative 
constants K and A, then for any constant p > 2, there exist a constant ci , depending 
on n and A such that 

P-2 

/ \ v \£*e- f dv] <- I (\Vy>\ 2 + r- 2 \p\ 2 )e- f dv 

JB (r) J V f (B (r))p JB (r) 

for allO <r < R and ipE C^(B (r)). 

Sketch proof of Lemma \Kl\ The proof is nearly same as that of Theorem 2.1 in [35] 
or Theorem 3.1 in |33j except our discussion with respect to the weighted measure 
e~'dv. When / is constant, this result was confirmed by L. Saloff-Coste [35] (see 
also Theorem 3.1 in [33 J. We refer the reader to these papers for the nice proof. □ 

Remark 2.5. In Lemma [2.41 the local Sobolev inequality is slight different from 
Munteanu and Wang's Neumann Sobolev inequality (Lemma 3.3 in [37]). Here 
we mainly follow the arguments of L. Saloff-Coste [32] to derive the local Sobolev 
inequality, whereas Munteanu and Wang proved their local Neumann Sobolev in- 
equality based on the argument in |16) . Also note that our local Sobolev inequality 
holds on any radius of balls due to a stronger assumption on /. But Munteanu and 
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Wang [27] established the local Neumann Sobolev inequality only on the unit balls 
due to a weaker hypothesis on the oscillation of / on unit balls. 

Below we shall present a result concerning the Harnack inequality for the /-heat 
equation, which is very much parallel to the case that / is constant, obtained by L. 
Saloff-Coste [32] and A. Grigor'yan [14]. 

Theorem 2.6. Let (M,g,e~ f dv) be an n- dimensional complete noncompact smooth 
metric measure space. Fix < R < oo. Assume that volume doubling property 
(|2.2p and the Neumann Poincare inequality (|2.3p are satisfied for this R. Then 
there exist constants C3 depending on n and A such that, for any ball B Q (r), o £ M , 
< r < R and for any smooth positive solution u of the f -heat equation in the 
cylinder Q = B a (r) x (s — r 2 , s), we have 



where Q- := B a (\r) x (s - |r 2 , s - \r 2 ) and Q + := B D (\r) x (s - \r 2 ,s). 

Sketch proof of Theorem \2.(A The proof is the same as the arguments of |32] or 
[33J except our discussion with respect to the weighted measure. Indeed this result 
follows from the standard Moser's technique. Since the conditions of Theorem 
12.61 imply a family of local Sobolev inequalities due to Lemma 12. 4[ combining the 
volume doubling property, this is enough to run the Moser's iteration procedure to 
prove Theorem 12.61 as explained in [32] or [33] . □ 

Combining Lemmas 12.11 12.21 and Theorem 12. 6[ we immediately have that: 

Corollary 2.7. Let (M, g, e~^dv) be an n-dimensional complete noncompact smooth 
metric measure space. If Ricf > —(n— 1)K and |/|(x) < A for some nonnega- 
tive constants K and A, then there exist a constant c(n, A) such that, for any 
< R < 00 and ball B Q (r), o £ M , < r < R and for any smooth positive solution 
u of the f -heat equation in the cylinder Q — B Q (r) x (s — r 2 , s), we have 



sup{w} < e c ^ A){1+Kr2) ■ inf{ U }, 
Q- Q+ 

where Q- := B Q {\r) x (s — |r 2 , s — \r 2 ) and Q + :— B (~r) x (s — \r 2 , s). 



In this section, we will apply the parabolic Harnack inequality to obtain a quan- 
titative Holder continuity estimate for a solution to the /-heat equation. Using 
this Holder continuity property, we shall derive a Liouville property under some 
suitable assumptions on Ricf and /. 

First, we give the Holder continuity estimate for any solution of the /-heat 
equation. When / is constant this was established in Theorem 5.4.7 of |34) . 

Theorem 3.1. Under the same assumptions of Theorem \2.6\ then there exist 9 = 



1 _ e -c(n,A)(i+Kr 2 ) g ^ ^ a = log2 g g ^ ^ and Ak = 4 0-l(l - > 1, where 

k G (0, 1) such that any solution u of the f -heat equation in Q = B (r) x (s — r 2 , s), 



sup{w} < e C3( - 1+Kr ^ ■ inf{u}, 



Q- Q+ 



3. Liouville theorem 




where Q 



K 



B {nr) x (s — nr 2 ,s). 
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Proof. The proof is nearly the same as the argument in [26] (see also [34]) which 
uses the parabolic Harnack inequality. For the reader's convenience, we include a 
detailed proof of this result. For any non-negative solution v of the /-heat equation 
in Q, by Theorem 12.61 we have 

(3.1) - 1 / vdp < max{u} < e < n ' A ^ 1+Kr ^ mhx{u}, 
Vf[Q-) Jq_ Q- Q+ 

where Q_ := B Q (\r) x (s — fr 2 ,s — \r 2 ) and Q + := B (^r) x (s — \r 2 ,s), and 
where Vf(Q—) denotes the volume of Q_ with respect to the space-time volume 
form dp. Now we let it be a solution, which is not necessary non-negative, and let 
M u , m u be the maximum and minimum of u in Q. Similarly, let M+, m+ be the 
maximum and minimum of u in Q + . Define 

where dp denotes the natural product measure on R x M: dp = dtx dp, and where 
dp = e~*dv. Applying (13.11) to the non- negative solutions M u — u, u — m u yields 

M u -^-<e< n ^ 1+Kr2 HM u -M+) 

and 

U, _ -m,< e« n ' A ^ 1+Kr2 \m+ - m u ), 

which imply that 

(M u - m u ) < e< n ^ 1+Kr2 \M u - m u ) - e «n,A)(l+Kr°) {M + _ m +y 
If we define the oscillations 

u)(u, Q) := M u - m u and u)(u, Q + ) := M+ - m+ 
of u over Q and Q + , then 

(3.2) w(«,Q+) < 9u(u,Q), 

where we assume e c(n,A)(i+Kr 2 ) > ^ and hence q = 1 _ e -c{n,A)(i+Kr 2 ) g ^ l). 
Now we consider (y,t), (y',t') £ Q K . Let 



p = 2 max{c%, y'), \Jt-t' } 

with t > t' . Then (y',t') belongs to Q := B y (p) x (t — p 2 ,t). We also define 
Pi = 2pj_x, po — p and Qi := B y {pi) x (t — pf,t) for all i>l. We easily see that 

(Qi)+ = Qi-l- 

Hence, as long as Qi is contained in Q, (13. 2[) yields 

uj(u,Qi-i) < 9uj(u,Qi) and u(u, Qo) < 9 l u(u, Q). 
Below, we consider two cases. If p < (1 — n)r, let k be the integer such that 

2 k < (l-«)r/p< 2 k+1 . 
Since (y,t) € Q K , it follows that 

Q k =B y (2 k p) x (i-4V,t) 

C - «)r) x (t - (1 - k) V, i) 

CB (r) x (s-r 2 ,s)=Q. 
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Hence we have 

w(u,Qo) < O k u(u,Q) < e-^l-K)-* (^j a u{u,Q) 
with a — log 2 0. This implies 

\u(y,t)-u{y',t')\ <^ sup{|u|} 



[{t-t'l 1 / 2 + d(y,y')] a ~ r a 

and conclusion follows, where A K := 46*~ 1 (1 — Kj~ a . The second case is trivial. 
Indeed if p > (1 — «)r, then the above inequality obviously holds. Therefore we 
complete the proof of theorem. □ 

Using the Harnack inequality and the above Holder continuity, we immediately 
derive following Liouville theorem. 

Corollary 3.2. Let (M, g, e~?dv) be an n- dimensional complete noncompact smooth 
metric measure space. Assume that Ricf > and \f\{x) < A for some nonnegative 
constant A. Then M has the strong Liouville property. Moreover, there exists an 
a € (0, 1] such that any f -harmonic function u which satisfies 

lim f r~ a ■ sup {\u\} I = 

for some fixed o G M is constant. 

Here, we say that a smooth metric measure space has the strong Liouville prop- 
erty if any solution u of the /-harmonic equation AfU — on (M, g, e~f dv) which 
is bounded below (or above) is constant. 



Remark 3.3. Corollary 13.21 was also proved by Munteanu and Wang 27 . We 
emphasize that this result can be regarded as a direct consequence of Theorem ll.il 
If / is constant and a — 1, then Corollary 13 . 21 returns to Cheng's Liouville property 
in [10]. If / is constant, this case appeared in [34] (see also Theorem 4.3 in [32 ). 



Proof of Corollarv \3. e A We first prove the first part of theorem. The conditions of 
theorem make the parabolic Harnack inequlity (Corollary 12. 7\ K = 0) hold and 
hence the corresponding elliptic Harnack inequality is also true. Assume that u is 
a solution of the /-harmonic equation which is bounded below. Let 

m(u) := infjuj. 

M 

Applying the elliptic Harnack inequality in the ball 2B = B Q (2r) and to the non- 
negative function v = u — m(u), we have that 

sup{u — m(u)} < C(n, A) ■ inf {it — m(u)}. 
b B 

We let the radius of B = B a (r) tend to infinity, and obviously infs{u — m(u)} 
tends to zero. Therefore we conclude that u — m(u) is constant. 

In the following we will prove the second part of theorem. Let a be as given by 
Theorem (5TJJ Let u be a function satisfying AfU — and 

lim ( r~ a ■ sup {M} | =0. 

r ^°° V B o(r) J 

Note that in this case we can choose a = 1. Because from the proof of Theorem 
13.21 it turns out that if 9 is sufficient small, then a — > 1. Fix some x £ M and y 
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such that d(x,y) < 1. Applying Theorem 13. II to u in a ball Br — B (R) with R so 
large that x, y € \Br, we find that 

C 

\u{x) - u(y)\ < -^sup{|u|}. 

Since the above inequality holds for all R large enough, we can let R tend to infinity 
to obtain that \u(x) — u(y)\ = 0. Since i,i/eM with d(x, y) < 1 are arbitrary and 
M is connected, we conclude that it must be constant. □ 

4. Two-sided Gaussian bounds on /-heat kernel 

In this section, we shall obtain upper and lower bound estimates for the /-heat 
kernel H(x,y,t) on complete noncompact metric measure space. The proof seems 
to be different from the classical discussion of Li and Yau in |21j . Our argument is 
similar to the discussion in A. Grigor'yan [14] and L. Saloff-Coste |32| . 

First, we show that the Neumann Poincare inequality and the volume doubling 
property imply a lower bound on the /-heat kernel. To achieve this, we start it by 
the following important lemma. 



Lemma 4.1. Under the same assumptions of Theorem \2.b\ then there exist a 
constant C4 := 04(71, A) such that, for any x,y £ B a (^r), and any < s < t < 00 
and any non-negative solution of the f -heat equation in M x (0, 00), 



1 1\ , . d 2 (x,y) 

K + — + -)(t-s)+ ; ,y) 



( u(x,s) \ 

Sketch proof of Lemma \4.1\ Since Theorem 12.61 implies a parabolic Harnack in- 
equality of the /-heat equation, it is sufficient to prove the above inequalities by 
carefully choosing different space-time solutions. Please see Corollary 5.4.4 in [34] 
or Corollary 5.4 in [33 for a detailed proof □ 

Using Lemma |4. 11 we can get a lower bound on the /-heat kernel on complete 
metric measure spaces. 



Proposition 4.2. Under the same assumptions of Theorem \2.b\ then there exists 
a constant C5 := 05(71, A) such that, for any x, y G B (hr) and any < t < 00, the 
f -heat kernel H(x,y,t) satisfies 



(4.1) H(x, y, t) > H(x, x, t) exp 



c 5 \l + — + Kt 



t 



Moreover, there exists a constant cq :— Ce(n, A) such that, for any x, y € B (^r) 
and any < t < r 2 < R 2 

-c 6 (l+Kt) / Mr \ 

(4.2) H(x, y, t) > exp -c 5 lX ' V) 



V f (B x (Vi)) "V t 

Proof. The proof of this result follows from that of Theorem 5.4.11 in [34] with 
no much modification. In fact using Lemma 14. 1[ we let u(y,t) = H(x,y,t) with 
x fixed and s = t/2 and then get (|4.ip . where we used the fact that H(x,x,t) is 
non-increasing. 

Then we prove (|4.2[) . Note that the conditions of this proposition imply a para- 
bolic Harnack inequality, which leads to the on-diagonal /-heat kernel lower bound 

(4.3) H(x,x,t) > e^^+^-Vf^B^Vt)) 
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for all x € M and < t < r 2 < R 2 . Indeed we fix < t < r 2 < R and consider <j> 
be a smooth function such that < <j> < 1, <j> = 1 on B := B x {yi) and <fi — on 
M\2B. Define 

' P t <t>{y) if t>0 

Hv) if <<o, 



u(y,t) 



Pt — e tAf be the heat semigroup of Af on L 2 {M,jj). Obviously, u(y,t) satisfies 
(dt — A/)it = on B X (—00,00). Using the parabolic Harnack inequality, first to 
u, and then to the /-heat kernel (y, s) — > H{x, y, s), we have 

l = u(x,0) <e c(1+Kt) u(x,t/2) 

= e c(1+m) / H(x,y,t/2)4,(y)dfx(y) 

JB(x,Vi) 

<e« 1+K V f H(x,y,t/2)d»(y) 

JB(x,2Vt) 

< e 2 < 1+Kt ^V f (B x (2Vi))H{x,x,t) 

< e M^ Vf (B x (Vt))e 2n+4A * Cln ' A)Vm H(x,x,t), 

where in the last inequality we used (|2 . 2[) . This gives (|4.3[) as desired. Hence (|4.2[) 
then easily follows by (|4.ip and (|4.3p . □ 



Secondly, we can show that the Neumann Poincare inequality and the volume 
doubling property also imply an upper bound on the /-heat kernel. To achieve this, 
the following integral estimate is critical useful due to E. B. Davies [llj . 

Lemma 4.3 (E. B. Davies [IT]). Let (M,g,e- f dv) be an n- dimensional complete 
smooth metric measure space. Let Ai > be the bottom of the L 2 -spectrum of the 
f -Laplacian. Assume that B\ and B2 are bounded subsets of M. Then 



f f H(x,y,t)d f i(y)d t i(x)<e- Xlt V f (B 1 ) 1 / 2 V f (B 2 ) 1 ^exp(~ 

JB l JB 2 \ 

where d(Bi,B2) denotes the distance between the sets B\ and B 2 . 



d 2 {B u B 2 ) 
At 



We are now ready to give an upper bound on the fundamental solution of the 
/-heat equation. 

Proposition 4.4. Under the same assumptions of Theorem \2.6l then there exist 
constants C7 and cs such that, for any x,y G B a (^r) and < t < r 2 /4, the f-heat 
kernel H{x, y, t) satisfies 

e c *<- 1+K ^ ( d 2 {x,y) 

Proof. Fix a fixed y £ B a (r) and 6 > 0, applying Lemma l4TT1 to the position solution 
u(x, t) = H{x, y, t) by taking s — t and t = (1 + S)t, 



H(x, y, t) < H(x', y, (1 + S)t) ■ exp <^ c 4 



I 1\ . , d 2 {x,x') 



K + — + -\ St 



St 
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Integrating over x' £ B x (y/i) gives 



H{x, y, t) < exp 



(4.5) 



c 4 (K + r~ 2 )5t + 5 + - 



H(x',y,(l + 8)t)dfi(x'). 
Applying Lemma 14.11 and the same argument to the position solution 

u(y,t) = H(x',y,t)dn(x'), 
by taking s = (1 + 8)t and t = (1 + 25)*, we obtain 



H(x', y, (1 + S)t)dfi(x') < exp 



B x {y/t) 



Substituting this into ()4.5j) yields 



c 4 (K + r- A )St + S + 



V f -\B y {y/t)) 



H(x, y, t) < exp 



H(x',y',(l + 2S)t)d^x')dn(y'). 



2c 4 (K + r- 2 )8t + 5 + - 



V-'iB^VtWf^ByiVt)) 



By (v^t) Jb x (VI) 
Combining this with Lemma |4.3[ we have 



H(x',y', (1 + 2S)t)dfi(x')dfi(y'). 



H(x, y, t) < exp 



(4.6) 



2c 4 ( (K + r- 2 )St + S+- ) -\ x t 



exp 



d 2 (B x (Vt),B y (Vi)) 
4(1 + 25)t 



Notice that if d(x,y) < 2\/i, then d(B x (y/t),B y (\/i) = and hence 



d 2 (B x (Vt),B v (Vt)) 



< 1 



d 2 (x,y) 



4(1 + 25)* 4(1 + 25)*' 

and if d(x, y) > 2 V*, then d{B x {-\fi), B v (y/i) = d(x, y) — 2\/t hence 

d 2 (B x (Vt),B y (Vt)) _ (d(x,y)-2V~tf < d 2 (:r,y) , 1 



4(1 + 25)* 4(1 + 25)* 

Therefore in any case, (|4.6[) becomes 



4(1 + 25)* 25 



H(x, y, t) < exp 



(4.7) 



2c 4 (K + r- 2 )St + S + 



Ai* 



d 2 (x,y) 



x ^ 1 / 2 (5 :c (v / *-))y / - 1/2 (i?,(v^))cxp (" 4(1 + M)t 
Now we want to estimate [K + r~ 2 )5* + 5 + ± in (|4.7p . Let 

5 = min{e ! [{K + r" 2 )*] ~ 1/2 } . 

If [(.ftr + r- 2 )*] -172 < e, then 

(K + r- 2 )5* + 5 + ~ < 2 [(X + r- 2 )t] 1/2 + e. 
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lf[(K + r~ 2 )t] ' >e, then we have 

1 
6 



(K + r- 2 )5t + 5+\<[{K + r- 2 )t] ( + 

1/2 1 



< [{K + r- 2 )t} 1/ +e + 



Hence, in either case, the right hand side of (|4.7[) can be estimate by 

1 

e 



H(x,y,t) < exp 



(4.8) 



d 2 (x,y) 
4(l + 2e)< 



2c 4 ^2 [(K + r 2 )f] 1/2 + £4 

x V f - 1/2 (B x (Vt))V f - 1/2 (B y (Vt))exp 
Moreover the volume doubling property implies (see, e.g., Lemma 5.2.7 in |34j ) that 
V f (x,V~t) < C(n,A)exp fc(n,A)VKi-^^j V f (y,Vi) 

< C{n, A) exp (c(n, A, e)Kt + V f (y, Vt). 

Substituting this into (|4.8[) and using < t < r 2 /4, then the theorem follows. □ 

Combining Lemmas 12. 1[ 12.21 and Propositions 14.21 14.41 we immediately obtain 
two-sided /-heat kernel bounds on smooth metric measure spaces. 

Theorem 4.5. Let (M, g, e~*dv) be an n- dimensional complete noncompact smooth 
metric measure space. If Ricf > — (n — 1)K and |/|(x) < A on B Q (2r) for some 
nonnegative constants K and A, then there exist positive constants Ci, i = 5, 6, 7, 8, 
depending on n and A such that the f -heat kernel H(x,y,t) satisfies 



3 (i+ift) 



Vf(x,Vt) 



exp —c, 



d 2 (x,y) 



<H(x,y,t) < 



e c 8 (l+Kt) 

v f (x,Vi) 



exp — C7 



d 2 (x,y) 



for any x,y € B a (r/2) and < t < r 2 /4. 



5. Mean value inequality 

In this section, the main objective is to derive a mean value inequality on com- 
plete noncompact metric measure space, which is a natural generalization of the 
Li-Schoen's result in [3D]. First, we give the following Poincare inequality. 

Theorem 5.1. Let [M, g, e~^dv) be a complete noncompact smooth metric measure 
space. If Ricf > — [n— l)K and \f\(x) < A for some nonnegative constants K and 
A, then for any a > 1, there exists constants C3 and C4 depending only on a, n 
and A such that 

I \4\°d» < C 3 f—^—Ye^VKr) I m „ dfl 
Jb {r) \1 + vKR/ Jb (r) 

for any compactly supported function <j> on B (R). In particular, the first Dirichlet 
eigenvalue \i\ of f -Laplacian on B Q (R) satisfies 
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Sketch proof of Lemma \5.1\ The proof is exactly the same as that of Corollary 1.1 
proved by Li-Schoen [20j except that the classical Laplacian comparison is replaced 
by the generalized Laplacian comparison (see Lemma |2.1[) 



A f r(x) < (n - 1 + 4A)V K cothV Kr 

< n - 1 + AA + (n-l + AA)VK. 
r 

Besides this, all the integration calculations should be done with respect to the new 
measure [i. To save the length of paper, we omit details of the proof. □ 

We now proceed to derive the L 2 mean value inequality with respect to the 
weighted measure by Theorem 15. II 

Theorem 5.2. Let (M, g, e~-> dv) be a complete noncompact smooth metric measure 
space. Assume that Ricf > — (n — 1)K with \f\(x) < A for some nonnegative 
constants K and A. Let u be a nonnegative f -subharmonic function defined on 
B (R). There exist a constant C§, depending only on n and A such that for any 
re (0,1/2) 



sup u 2 < t- c ^+^R)V7\B {R)) [ u 2 d^. 

B {{X-t)R) Jb (R) 



Proof. Let ft, be a harmonic function on B ((l - 2 1t )R) obtained by the solving 
the Dirichlet boundary problem 

A f h = on B ((1-t/2)R), 

and 

h = u on dB ((l ~ t/2)R). 

Since u is nonnegative, by the maximum principle, function h is positive on the ball 
B ((l - 2- 1 t)R). Moreover, 

u < h on B ((l - t/2)R). 

Using Lemmas 12.11 12 . 21 and 12 .41 by the Moser iteration argument as in [27], we have 
the following elliptic Harnack inequality 

SUp h < e c(n,A)d + VKR.) inf ^ 
B ((1-t)R) B ((1-t)R) 

where c depends only on n and A. In particular, 
sup u 2 < sup h 2 

B ((1-t)R) B {(1-t)R) 

< e c(n,A)(l+VKR) inf h 2 
\°- L > B ({1-t)R) 

< e c ^ A ^ 1+ ^^V7\B ({\ - t)R)) [ h 2 d[i. 

Jb ((i-t)r) 

Below we will estimate the L 2 (/x)-norm of h in terms of the L 2 (/i)-norm of u. By 
the triangle inequality, we have 

/ h 2 d[i <2 (h- u) 2 dn + 2 u 2 d^ 

^ 2^ Jb o ((1-t)R) Jb o ((1-t)R) Jb ((1-t)R) 

<2 / (h - u) 2 dfi + 2 1 u 2 dfi. 

Jb ((1-t/2)R) Jb (R) 
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Since (h — u) vanishes on dB {{\ — t/2)R) we can apply Theorem [5J] to show that 
(h-u) 2 d^< CsR2 e c * 1+VSs > [ \V(h-u)\ 2 dv 



B {(l-r/2)R) (1 + ^/KR) 2 JB ((1-t/2)R) 



q d2 dil+y/KR) r 

72- / 2 (|V^| 2 + |Vu| 2 )^, 

1 + VKR) JBo«1-t/2)R) 



where we have used the triangle inequality again. Since the Dirichlet integral of h 
is least among all functions which coincide with h on the boundary, from above we 
conclude that 

r 4f, Ft 2 p Ci(l + VKR.) p 

(5.3) / {h~u) 2 dii < 3 — / \Vu\ 2 dn. 

Jb o ((1-t/2)R) f 1 + ^/X R ) Jb o {{X-t/2)R) 



Now we use the fact that u is /-subharmonic to estimate the Dirichlet integral of 
u in terms of the Z/ 2 -norm of u. We have for any <j) with compact support in B Q (R) 

< / (J^uAfud/j, 



B (R) 

/ <f) 2 \\7u\ 2 dfi-2 4>u{Vcj),Vu)dn 
Jb (r) Jb (R) 



<-/ (t) 2 \Vu\ 2 dfi + 2[ [ cj) 2 \Vu\ 2 d^\ If u 2 \V(/)\ 2 dA . 

JB {R) \Jb o (R) J \Jb o (R.) J 

Thus 

f cj) 2 \\7u\ 2 dfi < 4 f u 2 \\7cf>\ 2 dfi. 

Jb o {R) Jb o {R) 

We let <f>[r{x)) be a cut-off function given by a function of r(x) = r(o,x) alone, 
such that 4>{r) = 1 on B ((l — t/2)R), 4>(r) = on dB (R), and satisfying 



™ * 7r- 



Then the above inequality becomes 



/ |V U | 2 d M < 44 / 

J BMI-t I2)R) T K J E 



■1 



u dfj,. 

,((l-r/2)fl) T~ IX~ J Bo (R) 

Combining this with (j5~Tj) . (|5~2)l and ([53)1 yields 
(5.4) 

sup u 2 <C 62CC3T % + 2 V f -\B D ((1 r)R)) / u 2 d^ 

B a {{l-r)R) \ (1 + VKR) 2 I JB (R) 



<c 6 T- c ^ 1 +^) e C8 ( 1+vl?fi )y- 1 (i? ((i - t)r)) J 



u 2 d/i 

B a (R) 



for some new constants d = Cj(n, A), i = 6,7, 8. To finish the proof, we also need 
to estimate the /-volume of B a (R) in terms of the volume of B ((l — t)R). Recall 
the bound for Ayr 2 : 

A/r 2 < 2(n + 4A) + 2\fK{n - 1 + 4A)r, 
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and hence 



/ A f r 2 dn<2(n + 4A)V f {t)+2VK(n-l + 4:A) / rdfi, 

JB„(t) JB a {t) 

where Vf(t) = Vol f(B a (t)). By Green formula, since 

A f r 2 dn= —do = 2t A ° ' 

s (t) yaB„(t) or dt 

then 

tV^(t) < (n + 4A)V)(i) + V^(n - 1 + 4A)tV f (t). 
Hence the function r!«+4A) e -v/x(r.-i+4A)i^ is decreasing in f > 0. Therefore 

(5.5) ^ / - 1 (i? ((l - r)i?)) < Vf X {B {R)) ( J • e ^Hr(n-l+4A) ) 

where < t < 1/2. Combining this with (|5.4p completes the proof of theorem. □ 

In the following, we show that the L p mean value inequality for any p G (0, 2] is 
a formal consequence of that given in Theorem 15.21 



Theorem 5.3. Under the same assumption of Theorem \ 5.2l for any p G (0,2], 
there exists a constant c depending only on n, p, and A such that 

sup u v <T-< lWSKi V7\R) f u p dn 

for any r G (0, 1/2), where Vf\R) := V7 1 {B {R)). 
Proof. By Theorem for any 5 G (0, 1/2], 6 G [1/2, 1 - 6], we have 
sup u 2 <S- c ^ 1+ ^^V f - 1 ((6 + S)R) [ u 2 dii. 

B (9R) JB ((e+S)R) 

Since 8 + 5 > 1/2, this inequality implies 

sup u 2 <S- c ^ 1 +^Vr 1 (2- 1 R) [ u 2 d^. 

B (6R) JB ((6+S)R) 

We also note that 

r ( Y~ p/2 r 

/ u 2 d[i < sup u 2 I u p d/j,. 

JB {(e+8)R) \B {(6+S)R) J JB ((8+8)R) 

Hence 

sup u 2 <5- c ^ 1+ ^^V7 1 (2- 1 R) ( sup u 2 ] f u"d^. 

B„(9R) ' \B ((6+8)R.) J J B ((6+8)R) 

If we set 

M(9) := sup u 2 

B (8R) 

and 



N :=V7 1 (2- 1 R) [ u p dfi, 

JBJR) 



(R) 

we have shown _ 

M(6) < NS- Cb( - 1+V1?r ^M(0 + S) 1 ^' 2 
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for any 6 e (0, 1/2] and 9 G [1/2, 1 — 5]. Choosing 

00 = 1 - t and 6i = 9i-i + 2~V 
for i = 1, 2, 3, we have that 

M(0i_i) < 7V 1 2 4C ' 5 ( 1+v/1? ^M(0 4 ) A , 
where A = 1 - p/2 and ATj = jtf T -C 6 (i+VWR) _ iterating yields 

M(0 O ) < if |t> W 2ft(iWjffl)!t 1 tt w J((« j )^ 
for any j > 1. Letting j tend to infinity yields 

M(0 O ) < T ~ C ^+^ R ) 

where Cg depends only on n, p and A. By the definition of M(9q), we have 

sup u p <T' 2 "P c ^ 1+ ^^Vf 1 (2- 1 R) [ uPdfi. 
B ((1-t)B) Jb (r) 

Finally, by the relation (|5.5p . i.e., 

V f - 1 (2- 1 i?) < C{n,A)e c{1+ ^ RR) Vf 1 {R), 

the theorem follows. □ 

6. L p -LlOUVILLE THEOREM 

In this section, we will study various L p Liouville theorems on complete noncom- 
pact smooth metric measure spaces. In the course of proving our results, integra- 
tions by parts, cut-off function techniques and comparison theorems are often used, 
which are more or less similar to classical results with respect to Ricci curvature 
case in [IB], [3D] and 03]. 

6.1. The 1 < p < oo case. The proof of this case is nearly the same as the 
classical case for the Ricci curvature proved by Yau [33] with little modification, 
but is included for completeness. 

Theorem 6.1. Let (AI, g, e~^dv) be an n- dimensional complete smooth metric 
measure space. Suppose o € M is a fixed point such that the geodesic ball B (2R) 
centered at o of radius 2R satisfies B Q (2R) n dM — 0. Let u be a nonnegative 
f -subharmonic function defined on B D (2R). Then for any constant 1 < p < oo, 
there exists a constant C(p) > depending only on p such that 

f U P~ 2 \Vu\ 2 dll < f U P d^. 

Jb o (R) R Jb o (2R.) 

In particular, if M has no boundary, then there does not exist any nonconstant, 
nonnegative, L p (p) -integrable f -subharmonic function. This constant must be zero 
if M has infinite f -volume. 

Proof of Theorem I6.il The proof is very similar to the Yau's arguments [43) . Let 
<p(r(x)) be a cut-off function given by a function of r(x) = r{o,x) alone, whose 
definition is as follows: 

r 1 if r < R, 
^ = {0 i ir >2R, 



-\2- l R) j 
Je 



u p dfi 



B„(R) 



4/p 
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with 4>{r) > 0, and < (4>') 2 /(f> < CR 2 for some positive constant C. Obviously, 
we have 

< / tfuP^Afudfi 

Jb o (2R) 

= -(p-l) I <p 2 u p - 2 \\7u\ 2 dn-2 I (j)U p - l {Vct),Vu)dii. 

JB (2R) Jb o (2R) 

Note that 



2 I (fyu^iVfrVu) < - I cfu p ~ 2 \Vu\ 2 + I u p |V</>| 2 . 



' B a (2R) 2 Jb o (2R) P~1Jb o (2R) 

Hence we have 

k - 1 



f u p - 2 \Vu\ 2 dn<^—^- f <j> 2 u p - 2 \Vu\ 2 dn 

•lB p (R) 2 Jb o {2R) 

< f u p \V(/}\ 2 dfi 

P- 1 JB (2R) 

- L 2R) uPdfi - 

If M has no boundary, by taking R oo, the fact that u is L p (/i) implies 

/ u p " 2 |Vm| 2 ^ = 0, 

Jb o {R) 

which furthermore implies |Vu| = 0. Hence u must be constant. □ 

6.2. The < p < 1 case. For < p < 1, following the work of Li-Schoen [20], we 
have the following result. 

Theorem 6.2. Let (M, g, e~ f dv) be an n- dimensional complete noncompact smooth 
metric measure space. If 

Ric f > -S(n)(l + r{x)y 2 , 

where r(x) is the distance to a fixed point o £ M , and f is bounded, then any 
nonnegative L p -integrable (0 < p < 1) f-subharmonic function must be identically 
constant. Moreover, this constant must be zero if M has infinite f -volume. 



Proof of Theorem \6.SX Since the arguments leading to Theorem 15.31 are local, by 
choosing more or less r = 1/2, we have the following L p mean value inequality 



(6.1) sup u p < ■2 c ^+V^(^ R )Vf 1 (B (R)) / vFdfi 

B {R/2) Jb o {R) 

for nonnegative /-subharmonic functions u on B X (5R), where Ricf > — (n — 
1)^(2;, 5i?) with |/|(x) < A for some nonnegative constants K and A on B x (hR). 
Here constant c depends on n, p and A. In the following, we will utilize (|6.1[) to 
show that u must vanish at infinity if nonnegative function u is /-subharmonic on 
M with u £ L p (n) (0 < p < 1). In fact, by the volume comparison theorem men- 
tioned above, under the hypothesis on Ricf and /, M must be of /-infinite volume 
and u must be identically zero. 
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Let x € M and consider a minimal geodesic 7 joining o to x such that 7(0) = o 
and j(T) = x, where T = r(o,x). We then define a set of values {U g [O,T]} k =0 
satisfying 

i 

t = o, t t = i + p, U =2^0* -l-F, 

where j3 > 1 to be chosen later, and tk — 2 53j=o — 1 — f3 k is the largest such 
value with < T. We denote the points Xi = 7(ii) and they obviously satisfy 

r(xi,x l+ i) = ft + /3' l+1 , r(o,Xi)=ti and r(x k , x) < f3 k + /3 fe+1 . 

Moreover, the set of geodesic balls B Xi (j3 l ) cover 7QO, 2 Ylj=o & ~ •"■]) an< ^ they have 
disjoint interiors. We now claim that 

/ on+iA \ k 

(6.2) V,{B mh {f)) >C[ W + 2 ^ +iA _ pn+iA ) V f (B (l)) 

for a fixed 

/3> 2V^t >l 

The proof of this claim essentially follows the arguments of Li-Schoen [5D] . For the 
sake of completeness, we will outline the proof of this claim again. 

For each 1 < i < k, a relative comparison theorem (see (4.10) in |39j ) argument 
shows that 

V f (B Xi (ft)) > A [Vf(B x ^ + 2/3- 1 )) - V f (B Xi (ft))] 
> DiVfiB^ift- 1 )), 

where 

D, = ■'" 



r(0 i +2P'-i)y/K(x i ,F +2F-i) sinh "-l+4A ^ ' 

since Ricf > — (n — l)K(xi,ft + 2/3 l_1 ) and |/|(x) < A for some nonnegative 
constants K and A on B Xi (/3 l + 2f3 t ~ 1 ). Iterating this inequality, we conclude that 

k 

(6-3) V f {B Xk (p k )) > V f (B (l)) JjA- 

i=l 

Since r(o, x,) = 2 X)j=o — 1 — the curvature assumption implies that 

(i-2 
J=0 

for i > 2. Hence 



P i y/K{x i ,p + 2/3 i ~i)< 



< 



2 2E;-= /3 j 



^■w-iKi-rr 1 
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which can be made arbitrarily small for a fixed /3 by choosing S(n) to be sufficiently 
small. Hence Di has the following approximation 

((3 l ) n+4A 



D, 



(p* + 2fr- 1 ) n + 4A - (/3 l ) n + 4A 



(ft _|_ 2)™+ 4 -4 _ fjn+iA 

by simply approximating sinhi with t. Hence (|6.2j) follows by combining (|6.3p . 

In the following, we shall estimate Vf(B x ((3 k+1 )). We achieve it by two cases. 
Case 1: r(x,Xk) < P k {P — 1)- In this case, we see that 

B Xk (p k ) C B x (p k+1 ), 

and hence we know 

V f (B Xh (P k ))<V f (B x (P k +')). 
Combining this with (|6.2I) . we conclude that 

/ an+iA \ k 

Vf{B x (P k +'))>C ^ + 2) P +iA _ pn+iA } V f (B (l)). 

Case 2: r(x,x k ) > P k (P — 1). In this setting, we see that 

B Xk {(3 k ) C B x (r(x,x k ) + p k )\B x (r(x,x k ) ~ p k ). 
Using a relative comparison theorem, we have that 

V f {B x (p k )) > D [V f (B x (r(x,x k ) + /3 k ) - V f (B x (r{x,x k ) - /3 k )] 
>D-V f {B Xk {p k )), 

where 

D = Jo 



\r(x,x k )-P k )y/ K(x,r(x,x k )+{t k ) 



Argument as above, since 



(r(x, x k ) + f3 k )^K(x, r{x, x k ) + /3 k ) < (p k+1 + 2/3 k )^K(x, r{x, x k ) + f3 k ) 
can be made sufficiently small, we can approximate D by 

pn+iA 

D ~ (/? + 2) n + 4A ' 
Combining this with (|6.2I) yields 

V f (B x (p k+1 )) > - +4A ( — - 4A _ p n+iA J V f (B (l)) 
~ ^ \ (p + 2) n + 4A - P n + 4A J V ^ B °^' 



where C depends on n, A and (3. In any case, (|6.4j) is valid. 
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If we let x — > oo, the value k — > oo. Note that the choice of 8 ensures that 

on+iA 

> 1 

(8 + 2) n+iA — f3 n + 4A ' 

and hence the right hand side of (|6.4|l tends to infinity. Now by (|6.1[) and the 
assumption of theorem, we have 

vP{x) < CV f -\B (R)), 

where C also depends on the L p -norm of u. Using the value R = B k+1 in the above 
inequality, the right hand side of the above inequality vanishes, thus proving that 
u(x) — > as x — > oo and Theorem 16.21 follows by the maximum principle. □ 

6.3. The p = 1 case. In the end of this section, we will consider the case p = 1. 
The proof of this case is a little complex. We shall follow the arguments of Li [TH] 
(see also [19] and [22]) to derive the following result. 

Theorem 6.3. Let (M, g, e~^dv) be an n- dimensional complete noncompact smooth 
metric measure space. Assume that 

Ricf > -C(l + r 2 (») 

for some constant C > 0, and \f\(x) < A for some constant A. Then any nonneg- 
ative L 1 (fx) -integrable f-subharmonic function must be identically constant. More- 
over, this constant must be zero if M has infinite f -volume. 

Following the trick of Li [18] (see also [22 ) to prove Theorem 16. 3[ at first, we 
need the following integration by parts formula. 

Theorem 6.4. Let (M,g,e~ f dv) be an n-dimensional complete noncompact smooth 
metric measure space. Assume that 

Ric f > -C(l + r 2 {x)) 

for some constant C > 0, and \f\{x) < A for some constant A. Then for any 
nonnegative L ((i) -integrable f-subharmonic function g, 

A f y H (x,y,t)g(y)dLi(y)= / H(x,y,t)A f g(y)dLi(y). 

M JM 



Proof of Theorem \6.4\ Similar to the arguments of [22., applying the Green formula 
on B (R), we have 

f A f H(x,y,t)g(y)dfx(y)-[ H(x,y,t)A f g(y)d Ll (y) 
Jb {r) Jb (r) 

[ ^-H{x,y,t)g{y)dLi aiR {y)- [ H(x,y,t)^-g(y)dfJ, (TtR (y) 
JdBJR) or JdBjR) or 



< \VH\{x,y,t)g{y)dLi a>R {y)+ H{x,y,t)\Vg\{y)dLi a>R {y), 

JdB (R) JdB (R) 

where Li a , R denotes the weighted are- measure induced by /z on dB {R). In the 
following we shall prove that the above two boundary integrals vanish as R — > oo, 
which can be achieved by four steps. 
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Step 1. Under the curvature assumption of Theorem 16. 4[ by the mean value 
inequality (see Theorem 15. 3p . there exist a constant c = c(n,A) such that 

(6.5) sup g(x) < e< 1+R ^V f - 1 (2R) [ g{yW{y). 

B (R) Jb o (2R) 

Consider cf)(y) = <t>{r(y)) to be a nonnegative cut-off function such that < <fi < 1, 
\V<j)\ < V3 and 

f onB (R+l)\B (R), 
on B (R-l)U (M\B (R + 2)). 

Since g is /-subharmonic function, by the Schwarz inequality we have 

0< / c/) 2 gA f gdpi = - / V {(j) 2 g)V gd/i 
Jm Jm 



<j>(x) 



= -2 f 4>g{V(jNg)dn- [ </> 2 |V. 9 | 2 ^ 

<2 f |V0|VdAt- | ! 2 |V 5 | 2 dM. 
Jm 1 Jm 



im * Jm 

Then using the definition of <j) and (|6.5p . we have that 



|V 5 | 2 ^ <4 f \V^\ 2 g 2 d^ 
Jm 



B (R+1)\B (R) JM 



<12 / g^d/i 

>B (R+2) 



<12 sup g ■ / gd/i 

B (R+2) JB (R+2) 
12e c(l+RVK) f , \ 2 

*V f {2R + 4) ■{L R+ 2 ) 9dfl 



From this, using the Schwarz inequality, we get 

/ \ 1/2 

/ \Vg\dfx < j |V.g| 2 d M ■ [V f (R + l)\V f (R)}^ 2 

J B (R+1)\B (R) \J B (R+1)\B (R) J 

<( I |V 5 | 2 d^ ■V f (2R + 4) 1 / 2 . 

\Jb (R+1)\B (R) J 

Combining the above two inequalities, we have 

(6.6) / \Vg\dn < C 10 e< 1+R ^ [ gdfi, 

J B (R+1)\B (R) Jb o (R+2) 

where Cio = Cio(n, A). 

Step 2. We estimate the /-heat kernel H(x,y,t). By Theorem 14.51 the /-heat 
kernel H (x, y, t) satisfies 

e^ 1+Kt "> ( d 2 {x,y) 

(6-7) H(x,y,t) < — — -=- exp -c 7 



V f (B x (Vt)) "\ 1 t 
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for all x,y £ B {R) and < t < R 2 /8, where the constants C 7 and Cs depending 
on n and A. Now combining (|6.6p with (|6.7|) gives 



Ji ■■= 



B {R+X)\B {R) 



H(x,y,t)\Vg\(y)d»{y) 



< sup H(x,y,t) 

\yeB (R+i)\B (R) 



< — • exp 



'V f (B x (Vt)) 



V f {B x {Vi)) Jb (r+2) 



c s Kt- 
fJ 



B (R+1)\B (R) 

c 7 (R-r) 2 ^ 



t 



gdn 



B (R+2) 



x exp 



cVKR + c 8 Kt 



c 7 R 2 2c 7 Rr c 7 r 2 



t 



t 



t 



where r — r(o,x) and C\\ = Cn(n,A). We now confirm that J\ — > as R — >• 00. 
Indeed, for any fixed point x £ M, if we choose time t sufficiently small fixed such 
that — dominates a major part with respect to the quadratic inequality in R, 
then Ji tends to zero as R tends to infinity. 

Step 3. Below we shall estimate the gradient of H. We consider the integral 
with respect to dfi: 



M 



</» 2 (j/)| VH\ 2 {x, y, t) = -2 / (H(x, y, t)V0(y),0(y)Vff (i, y, t)) 



A I 



M 



2 (y)H(x,y,t)A f H(x,y,t) 



<2 \Vcj>\ 2 (y)H 2 (x,y,t) + - <f> 2 (y)\VH\ 2 (x,y,t) 



A I 



A I 



cf> 2 {y)H(x,y,t)A f H(x,y,t). 



This implies 



B (R+1)\B (R) 



< 



cf> 2 (y)\VH\ 2 (x,y,t) 



M 



< 4 / \V<t>\ 2 H 2 -2 (j> 2 HA f H 



< 12 



< 12 



M 

H 2 + 2 



B (R+1)\B (R-1) 



B (R+2)\B (R-1) 



H\A f H\ 



H 2 + 2 



,(R+2)\B (R-1) 



H 



B (R+2)\B (R-1) 



A I 



We point out that it is known that the heat semi-group is contractive in L , hence 

H{x,y,t)dfi{y) < 1. 



A I 
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Using this and (|6.7|) . we can estimate 



H (x,y,i)dfj, < sup H(x,y,t) 

J B {R+2)\B (R-1) ' yeB (R+2)\B (R-X) 

( 6 - 9 ) ^ r „ _r_ „ M 2 



< 



C\2 



Vf(B x (Vt)) 



x exp 



c 8 Kt 



-c 7 (R-r(o,x)y 



t 



Also, we claim that there exists a constant C13 > such that 
(6.10) 



f (A f H) 2 (x,y,t)d»<^-H(x,x,t). 



I M 

To prove this inequality, we first derive the inequality for any Dirichlet /-heat 
kernel H defined on a compact subdomain of M. Using the fact that /-heat kernel 
on M can be obtained by taking limits of Dirichlet /-heat kernels on a compact 
exhaustion of M, then (|6. 10[) follows. Indeed, if H(x,y,t) is a Dirichlet /-heat 
kernel on a compact subdomain fl C M, using the eigenfunction expansion, then 
H (x, y, t) can be written as the form 



i 

where {V'i} are orthonormal basis of the space of L 2 functions with Dirichlet bound- 
ary value satisfying the equation 

Afipi = -Xiipi. 

Differentiating with respect to the variable y, we have 

00 

A f H(x,y,t) = -^Ke-^^x^iy)- 

i 

Noticing that s 2 e~ 2s < C\^e~ s for all < s < oo, therefore 

p oo 

/ (AfHfdniy) < C 13 t- 2 ]T e-**'tf{x) 
Jm 

= C 13 r 2 H(x,x,t) 

and claim ([6~TU|) follows. Now combining (1B~H) . (|6~9)) and ([6~TU| . we obtain 
/ \S7H\ 2 dfi <C*i4 Wf 1 +t- 1 V7^H^{x 7 x 7 t) 

J B (R+1)\B (R) 



x exp 



c 8 Kt 



-c 7 {R~r(o,x)f 
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where Vf := Vf(B x (>/t)). Applying Schwarz inequality and volume comparison 
theorem, we have 



B (R+1)\B (R) 



\VH\dp < [V f (B (R+l))\V f (B (R))] 



1/2 



(6.11) 



1 V2 



\VH\ 2 dfi 



B (R+1)\B (R) 

-1 , 



<C(n,K,A) Vf +t~ Vf 2 H^(x,x,t) 



1/2 



x exp 



2 2t 



for sufficiently large R > 0. On the other hand, by (|6.5|) . (|6.11|) and Schwarz 
inequality we see that 



Ji ■= 



\VH(x,y,t)\g(y)d^(y) 



B (R+1)\B (R) 



< sup g(y) 

yeB (R+l)\B (R) 

C 15 e c ^ R 



\VH(x,y,t)\dti(y) 



Ba{R+l)\B {R.) 



< 



V f (B (2R + 2)) J B 
x exp 



C ±Kt 
2 



g(y)dfi(y)- Vf 1 +t~ 1 V f *H*(x,x,t) 

(2R+2) L 

-c 7 (R- d(o,x)) 2 



1/2 



2t 



Similar to the discussion of J%, for any a fixed point x € M, and any sufficiently 
small fixed t > 0, J 2 also tends to zero when R tends to infinity. 



Step 4- Recall that the co-area formula states that for all h 6 Cq°(M) 

r rR+l 

/ Ky)dfi(y) = 

Je 



I B (R+1)\B (R) 



dB (R) 



B o(r)h(y)dfj, a!r (y) 



dr. 



By the mean value theorem, for any R > there exists R € (R, R + 1) such that 



[H(x, y, t)\Vg\(y) + \VH\(x,y, t)g(y)] d^ R (y) 

8B (R) 

[H(x,y,t)\Vg\(y) + \VH\(x,y,t)g(y)] d^y) 



J : = 



B (R+1)\B P (R) 
= Jl + J 2 - 

By step 1 and step 2, we know that for any a fixed point x € M, and for any sufficient 
small fixed t, J tends to zero as R (and hence R) tends to infinity. Therefore we 
complete Theorem 16.41 □ 



Now we can finish the proof of Theorem 16.31 

Proof of Theorem \6.3[ Let g be a nonnegative, L 1 -integrable and /-subharmonic 
function defined on M. Now we define 



g(x,t) 



H (x,y,t)g(y)d^i(y) 



M 
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with g(x, 0) = g(x). By Theorem I6.4[ 
d f d 

Qj.9(x,t) = J —H(x,y,t)g(y)diJ,(y) 
A f v H (x,y,t)g{y)dn(y) 

M 

H(x,y,t)A f g(y)dn{y) > 0. 

M y 

Therefore g(x, t) is increasing in t. Meanwhile, by the mean value theorem, for any 
R > there exists R £ (R, R + 1) such that 



Since 



[ \VH\(x,y,t)d t i aM {y)= [ |V H\{x,y,t)dfi{y). 

JdB (R) J B (R+1)\B (R) 

[ A f H(x,y,t)dfjL(y) < f \VH\(x, y,t)d^ R {y), 

JBJK) JdBJR) 



so we have 



/:=/ A f H(x,y,t)dfi(y) 
Jb (r) 

< [ \VH\(x,y,t)dn(y) 

J B (R+1)\B (R) 

1/2 



<C(n,K, A) 



Vf 1 +t- L V f 2 Hi(x,x,t) 



x exp 



—Kt < ~ c 7( i? " d (°> a; )) 2 



2 2t 

for sufficiently large R, where we used (|6.1ip for second inequality above. Similar 
the above arguments, for any a fixed point x £ M, if we choose t sufficiently small, 
then / tends to zero when R and hence R tends to infinity. Hence 

| - [ H(x,y,t)dfi(y)= [ A f H(x,y,t)dfi(y)=0, 
at Jm Jm 

which implies 

H(x,y,t)dfi(y) = 1 



IM 

for all x £ M and for all sufficient small fixed t. Moreover, this equality implies 
g(x,t)dx= / H(x,y,t)g(y)dii(y)dii(x) = g{y)dy. 

M JMJM JM 

Since g{x,t) is increasing in t, we conclude that g(x, t) = g(x) and 
A f g(x) = A f g(x, t) = [ A f H(x, y, t)g{y)dfi(y) 

= 5(0- / A f H{x,y,t)dfx(y) = 0, 

JM 

On the other hand, for any positive constant a, we define a new function 

h(x) :— min{g(a;), a}. 

Then it satisfies 

h{x) < g(x), \Vh\ < \Vg\ and A f h(x) < 0. 
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In particular, it will satisfy the same estimates, (|6.5[) and (|6.6|) . as g. Hence we can 
show that 



Note that h is still L 1 , following the same argument as before, we have Afh(x) = 0. 
Hence regularity theory asserts that h must be smooth. Since a is arbitrary, this is 
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